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Abstract Understanding the subgap behavior of Normal-Insulator-Superconductor 
(NIS) tunnel junctions is important in order to be able to accurately model the 
thermal properties of the junctions. Hekking and Nazarov ' developed a theory in 
which NIS subgap current in thin-film structures can be modeled by multiple An- 
dreev reflections. In their theory, the current due to Andreev reflections depends 
on the junction area and the junction resistance area product. We have measured 
the current due to Andreev reflections in NIS tunnel junctions for various junction 
sizes and junction resistance area products and found that the multiple reflection 
theory is in agreement with our data. 
PACS numbers: 74.45 .-i-c 74.50.-i-r 



1 Introduction 

Accurate modeling of the current-voltage (IV) characteristics of NIS junctions is 
required in order to use the junctions in applications such as primary thermometers 
or solid-state refrigerators. The Bardeen-Cooper-Schrieffer (BCS) theory of super- 
conductivity accurately predicts the IV characteristics of NIS junctions above the 
superconducting energy gap A, while predicting almost no current when the junc- 
tion is biased below the gap^ ^. When NIS junctions are measured, subgap currents 
greater than BCS predictions are often measured. This excess current can be ex- 
plained by Andreev reflections^^, where an electron (hole) in the normal metal is 
reflected from the NS interface as a hole (electon), which allows a Cooper pair to 
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enter (leave) the superconductor. In interfaces where the electrons and quasiparti- 
cles are in the ballistic regime and can be represented as a plane wave, the Andreev 
reflection is described by the Blonder- Tinkham-Klapwijk (BTK) theory^. This ex- 
cess power load in NIS junctions was modeled by Bardas and Averin^. In realistic 
interfaces, the electrons and quasiparticles are no longer in the ballistic regime 
but behave diffusively, because they can reflect off the barrier and surrounding 
surfaces many times before tunneling, which will cause a higher current than pre- 
dicted by the BTK theory. Hekking and Nazarov developed a model to account for 
the extra current due to multiple Andreev reflections. Rajauria et al. ^ have made 
measurements of SINIS IV curves, and their subgap data agreed with Hekking 
and Nazarov's theory when they multiplied their data by a scaling factor of 1.371^. 
However, Rajauria et al. preformed measurements on only a single junction size 
and a single oxidation thickness. In order to understand NIS junctions, it must be 
understood how junction area and interface resistance affect the Andreev current. 
In this paper, we provide a more robust test of Hekking and Nazarov's theory by 
comparing it with measurements of multiple NIS devices with different junction 
areas and resistance area products. 



2 Theory 

BCS theory predicts that current will flow through an NIS junction at voltage bias 
greater than the superconducting energy gap. When the bias voltage is less than 
the superconducting gap, little current should flow because electron tunneling is 
limited by the density of states in the superconductor. However, a current below 
the gap is possible from mechanisms such as Andreev reflections. Hekking and 
Nazarov predict an additional current below the subgap caused by Andreev reflec- 
tions: 



hndreev ~ ^5"^ tanh(eV'/2r) + -5-^^ , (1) 

e^RlSVNdN ^ ' ' e^RlSVsds2KA^\-eV/A 

where e is the electron charge, /?„ is the normal state resistance of the junction, S 
is the junction area, Vns is the density of states of the normal metal (superconduc- 
tor), c/yv.s is the thickness of the normal metal (superconductor), V is the voltage 
bias and T is the temperature. By rearranging Eq.[T| we find that the dimensionless 
quantity eIR„/A scales inversely with the resistance area product. 

Equation [T] is valid as long as the junction dimensions are larger than the 
coherence length in each metal. The coherence length of the superconductor is 
given by ^5 = y^hDs/A, and the coherence length of the normal metal is given 

by = y/'hDN /k},T , where Ds,n is the diffusion constant for the superconductor 
and the normal metal, respectively. We measured the resistivity of our supercon- 
ductor, Al, to be Pai = 0.005 Q x fim, and measured the resistivity of our normal 
metal, AlMn'", to be pAiMn = 0.10 Q x /im. Using these values, we calculated 
the coherence lengths to be ^5 = 444 nm and = 511 nm. Our smallest junction 
dimension is 2 fim, which is greater than both coherence lengths. Therefore, it is 
valid to use Eq.[T]to model our junctions. 
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Fig. 1 (Color online) Optical image of a device used in our measurements. The junction is made 
by overlapping normal metal and superconducting wires. The junction area is defined by a via 
in a Si02 layer between the two metal layers. 

3 Experimental Details 

Figure [T] shows a typical device used in this experiment. The junctions were cre- 
ated by intersecting normal metal and superconducting wires, and the junction 
area is defined by a via in a layer of Si02 separating the two metal layers. The 
junctions were fabricated on a Si wafer by first sputter depositing 30 nm of Al 
doped with Mn to 4000 ppm by atomic percent. The AlMn was patterned using 
standard photolithographic techniques and was etched in an acid bath. A 90 nm 
thick layer of Si02 was deposited by use of plasma-enhanced chemical vapor de- 
position (PECVD) and vias were created by use of a plasma etch to define the junc- 
tion area. The devices were ion milled to remove any native oxide from the normal 
metal, and then exposed to oxygen to form the insulating layer. A resistance area 
product, Rsp = R„ x S, of 30 X2/im^ was created by exposing the devices to 0.1 
torr-s of oxygen and a Rsp of 200 12 /im^ was created by exposing the devices to 
42 torr-s of oxygen. Finally, the superconducting Al counter electrode was sputter 
deposited and then wet etched by use of standard photolithographic techniques. 
The 30 12 /im^ devices were fabricated on top of a 150 nm layer of Si02 deposited 
by PECVD to increase the quality of the junctions, but no observable difference 
was measured between these junctions and junctions that were fabricated on just 
the thermal oxide. 

The junctions were then screened to determine their quality before measure- 
ments were performed. To determine the quality of our junctions, we use the qual- 
ity factor Q, where Q = Rieak/Rn- The leakage resistance, Rieak, is defined as the 
highest resistance of the junction in the subgap. Both wafers produced devices 
with a e w 2000. 

For the experiment, we measured four devices, two from each wafer we fab- 
ricated. We measured devices 3 /im by 3 /im and 4 /im by 4 /im from the wafer 
with a Rgp = 200 12 /im^, and devices 2 ^m by 2 /im and 3 /im by 3 /im from the 
wafer with a Rsp = 30 12 /im^. Devices from the same wafer were chosen from the 
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Device 


Area (/xm'^) 


R„ (f2) 


Rsp (n^m^) 


djv (nm) 


ds (nm) 


A (/xeV) 


1 


16 


11.2 


179 


30 


525 


185 


2 


9 


21.5 


193 


30 


525 


185 


3 


9 


2.98 


27 


31 


230 


185 


4 


4 


8.25 


33 


31 


230 


185 



Table 1 Measured values of device parameters used in our experiment. 
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Fig. 2 (Color online) Current vs voltage for the four devices that we measured. The data are 
represented by points and the theory as dashed lines. The error bars for the uncertainty in the 
measurement are smaller than the data points. The theory was calculated using the device pa- 
rameters shown in Table 1 . The data are in excellent agreement with the theory at lower voltage 
biases. At higher biases, a detailed thermal model is needed to match the theory with data. 



same chip to make sure that the junction properties, such as the metal thicknesses 
and Rsp, were as similar as possible. 

IV measurements of the devices were made by four wire measurements in an 
adiabatic demagnetization refrigerator at 100 mK. Current biasing was accom- 
plished with a low-noise voltage source in series with a 10 MQ resistor. Table 1 
shows the properties of the devices that were measured in this experiment. The 
normal state resistance was measured from the differential resistance of the de- 
vices. The normal metal and superconductor thicknesses were measured for each 
device by use of a profilometer. 



4 Results 

The results of our experiment are shown in Fig.|2] The uncertainty in the measure- 
ment is ± 50 pA and ± 5 /i V, which makes the error bars smaller than the markers 
in the plot. The current was scaled by eR„/A, and the voltage was scaled hy e/A 
in order to make the axes unitless. The theory lines were calculated by using Eq. 
[T| with the measured device parameters shown in Table 1 . The divergence of the 
data with the theory at higher voltage biases, around 0.7 A je and above, is due to 
the non-isothermal behavior of the NIS junctions. In this work, we modeled the 
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Fig. 3 (Color online) Current vs voltage for the first two devices plotted against NTS curves 
generated using various Dynes parameters. The Dynes theory does not fit the data in the subgap. 
The Dynes theory predicts a steadily rising current in the subgap, which we do not see in our 
data. Instead, we see a flat plateau below the subgap, which is consistent with Andreev reflection. 



devices as being isothermal. However, NIS junctions are known for their ability 
to cool electrons in the normal metal at bias voltages near the superconducting 
gap^. Therefore, in order to calculate the true behavior of these junctions, the IV 
curves need to be modeled by solving a complex power balance equation. Since 
we are interested only in current due to Andreev reflections, which occurs below 
the gap, the junctions behave isothermally in the region of interest and no power 
balance modeling is required. 

As Fig. |2] shows, our data are in excellent agreement with the theory below 
the superconducting gap. Devices with the same Rsp have the same dimensionless 
current, and the current due to Andreev reflections is inversely proportional to Rgp, 
as theory predicts. Dividing the subgap data by the theory. Devices 1 and 2 agree 
with theory within 7 %, Device 3 agrees within 1 8 % and Device 4 agrees within 
15%. 

For comparison, we also fit our data with NIS IV curves based on the phe- 
nomenological Dyne's parameter'^. The Dynes parameter is added to the BCS 
density of states in an attempt to account for the broadening of the gap edge and/or 
the presence of subgap states, as shown in Eq. |2] 
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Re 



E-ir 
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dE (2) 



where /n is the Fermi function of the normal metal and F is the Dynes parameter. 

As Fig.|3]shows, the Dynes theory does not provide a good fit to our data. The 
Dynes theory predicts a steadily rising current below the gap, while in our data, 
the current plateaus, which is consistant with Hekking and Nazarov's theory. This 
supports that we are measuring Andreev reflections and not subgap conductance 
due to the presence of subgap states. 
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5 Conclusion 

In this paper, we have measured the current due to mukiple Andreev reflections in 
NIS junctions and found that it is in excellent agreement with the theory presented 
by Hekking and Nazarov. Their theory is valid only for voltage biases below the 
superconducting gap, and we are working to incorporate their theory with our 
power balance equations^^lin order to make more accurate thermal models of NIS 
junctions. These models will allow us to better predict the cooling properties of 
junctions with a small Rsp, for which the Andreev current is significant. 
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